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Abstract 

Based on a microscopic density functional theory we calculate the internal 
structure of the three-phase contact line between liquid, vapor, and a con- 
fining wall as well as the morphology of liquid wetting films on a substrate 
exhibiting a chemical step. We present a refined numerical analysis of the non- 
local density functional which describes the interface morphologies and the 
corresponding line tensions. These results are compared with those predicted 
by a more simple phenomenological interface displacement model. Except for 
the case that the interface exhibits large curvatures, we find that the inter- 
face displacement model provides a quantitatively reliable description of the 
interfacial structures. 
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I. INTRODUCTION 



The formation of a homogeneous wetting film of phase (3 with thickness Iq at the planar a- 
7 interface between two coexisting phases a and 7 is governed by the interplay of the surface 
tensions cr a p, <7g 7 , and a ai 0^1] • In the standard case (3 and 7 are the liquid and vapor 
phases, respectively, of a simple fluid and a represents a substrate acting as an inert spectator 
phase. The translational invariance of these systems in the lateral directions can be broken 
either by opposing boundary conditions or by geometric or chemical heterogeneities within 
the confining medium a. Whereas the first mechanism can always be applied, irrespective 
of the nature of the phase a, the latter two mechanisms require a solid phase a which can 
permanently sustain well-defined lateral structures. 

As a prerequisite for the possibility to impose opposing boundary conditions in lateral 
directions the thermodynamic state of the system has to be chosen such that the two phases 
/3 and 7 are in thermal equilibrium. In that case the system can be arranged such that on 
its left end the substrate a is exposed to the vapor phase in the bulk whereas on its right 
end the substrate is exposed to the coexisting liquid phase in the bulk. This arrangement 
leads to the formation of a liquid- vapor interface intersecting the substrate at a three-phase 
contact line with the contact angle 6 = arccos((cr^ 9 — a w i)/crig) (see Fig. |l|). The free energy 
of this configuration decomposes into the volume contributions corresponding to the gas 
and liquid phases, respectively, the surface tensions a wg , a w i, and a\ g associated with the 
corresponding half-planes of contact between these phases, and the line tension r associated 
with the presence of the three-phase contact line |I],|5|,f|. 

There have been numerous theoretical (see, e.g., Refs. |7Hl5f) and experimental (for a 
review see, e.g., Ref. |]16[) efforts to determine the sign, magnitude, and temperature depen- 
dence of line tensions (see Refs. |l7j and |18|] for a recent summary of an extended list of 
additional references). In particular the singular behavior of the line tension upon approach- 
ing wetting transitions, i.e., for 9 — > has been examined by using simple phenomenological 



interface displacement or square-gradient models [p^,|l9|-|2"2|j . Since the long-ranged, i.e., 
power-law decay of the dispersion forces acting in fluids is known to invalidate a gradient 
expansion for effective interface Hamiltonians [T2^-pB[], in Ref. |I7| the actual nonlocal inter- 
face Hamiltonian, as it is obtained from a microscopic nonlocal density functional theory, 
has been used to study the influence of dispersion forces on the line tension and on the 
intrinsic structure of the three-phase contact line. This study led to significant quantitative 
differences in comparison with the more simple interface displacement model. 

Similar differences between the predictions of the nonlocal theory and its local approxi- 
mation appeared []27] in the analysis of the morphology of a wetting film covering a planar 
substrate which consists of two adjacent halves composed of different materials and thus 



represents a chemical step (see Fig. g). According to the results reported in Ref. |27] the 
nonlocal theory predicts a much broader transition region, within which the local film thick- 
ness l(x) switches as function of the lateral coordinate x between its asymptotic values 
l± = l(x — > ±00), than the interface displacement model does. 

The local interface displacement and square gradient theories relish popularity and are 
convenient for the description of thin fluid films on chemically or geometrically structured 
substrates (see, e.g., Refs. [p8|-|39|j). In views of the aforementioned known limitations of the 
applicability of square gradient theories for systems governed by dispersion forces |23|-p6| the 
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objective of the present study is to pinpoint the reason for the reported large quantitative 
differences following from the local and the nonlocal approach. The common expectation is 
that the interface displacement model and the square gradient theory, although they are only 
approximations to the full nonlocal theory, turn out to yield nonetheless reliable results for 
most of the cases. The studies of the three-phase contact line (Sec. |IT|) and of the adsorption 
on a chemical step (Sec. [HID serve as testing grounds for the comparison between the local 
and the nonlocal theories. Our results are summarized in Sec. ffVI. 



II. THREE-PHASE CONTACT LINE 



As stated in the Introduction we consider a simple fluid in a grand canonical ensemble 
described by the chemical potential p and the temperature T. p = po(T) is chosen such 
that in the bulk the fluid is at liquid-vapor coexistence. This allows one to maintain the 
configuration shown in Fig. |I] which is characterized by the number density p(x, z) of the fluid 
particles. It interpolates smoothly between p(x — > —00,2) = p~(z), which corresponds to 
the wall- vapor interfacial profile, and p(x — > 00, z) = p+(z), which describes the wall-liquid 
interfacial structure. The local position l(x) of the liquid- vapor interface can be determined, 
e.g., as the isodensity contour line p(x, z = l(x)) = \{pi + p g ), where pi and p g are the bulk 
densities in the liquid and vapor phase, respectively. l(x) asymptotically approaches the 
finite value Iq(T) for x — > —00, which corresponds to the equilibrium wetting film thickness 
at the wall-vapor interface, and diverges linearly for x —>■ 00 with a slope given by the 
contact angle (see Fig. |l|). 



A. Density functional theory 

Density functional theory has turned out to be a suitable theoretical description for 
spatially inhomogeneous fluids as considered here. We apply a simple [3D] but nonetheless 



successful version which captures the essentials of wetting transitions ||. Its grand canonical 
free energy functional reads 

n([p(r)};T,p)= [ d 3 rf HS (p(r),T) + [ d 3 r[V(r) - p}p(r) 

J A J A 

+ - / d 3 r [ d 3 r'w{\r-r'\)p{r)p(r'). (2.1) 

2 J A J A 

A is the finite volume filled by fluid within the half space A + = {r e M. 3 \z > 0}. In the 
thermodynamic limit one has A — > A + . The other half space {r e M. 3 \z < 0} is occupied by 
the substrate. The external potential V(r) describes the interaction of a fluid particle with 
the substrate: 



V(r) = V(z>0) = -J2^- (2-2) 



i>3 



Approximately V(z) can be thought of as a laterally averaged pairwise sum of Lennard- Jones 
potentials 
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between individual fluid and substrate particles. 

w(r) describes the attractive part of the pair potential between the fluid particles which 
is taken to be of Lennard- Jones type, i.e., <fif(r) = 4e/[(<T//r) 12 — (cr//r) 6 ]. A division scheme 
like the Weeks-Chandler-Andersen (WCA) procedure provides an attractive {4>att) an d 
a repulsive part (0 rep ) as suitable entries into the expression (|2.1|). We approximate the 
attractive part by the smooth function 
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(2.5) 



w(r) exhibits the main feature of the attractive van der Waals interaction, namely the large- 



distance behavior w(r) 



Its amplitude is chosen such that the integrated strength of 



w(r) is the same as the integrated strength of 4>att as constructed from the WCA procedure. 
The repulsive part 4> rep (r) = 6(2 1 / 6 o"/ — r)((f)f(r) + £/) (with the Heaviside step function 0) 
of the pair interaction gives rise to a reference free energy fns{p,T) of a hard sphere fluid, 
for which we adopt the Carnahan-Starling approximation f4"2fl: 



f HS (p,T) = k B Tp ln( P A 3 )-l 



4r] — 3?7 S 



(2.6) 



with the thermal de Broglie wavelength A, the dimensionless packing fraction rj = |p(c/(T)) 3 , 
and the effective hard sphere diameter 



2 1 /6 (T , 



d(T) 



dr < 1 — exp 



■'rep \ 
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(2.7) 



In Eq. ( |2.1| ) the reference free energy is evaluated in a local density approximation and 
therefore it does not properly take into account the details of the packing effects near the 
wall. If one would be interested in these aspects, more sophisticated density functional 
theories have to be applied. Due to the last term in Eq. ( [2.1|) the present density functional 
is a nonlocal expression. The bulk phase diagram, i.e., the values for the bulk particle 
densities pi(T, /i) and p g (T,p) can be calculated by minimizing the bulk free energy density 



Q b (p,T,p) 



fiis(p,T) + -w p 2 



pp 



(2.8) 



with respect to p, where wo is defined in Eq. ( |2.5| ). Equation (|2.8| ) follows from inserting the 
homogeneous bulk density into Eq. (|2.1|). At two-phase coexistence one has 



dp 



dn h 



p=pg 



dp 



and fl b (Pg) = &b{pi 



(2.9) 



p=pi 
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With the bulk properties fixed the functional expression (|2.1|) can now be used to analyze 
the morphology and the line tension of a three-phase contact line. In spite of the relative 
simplicity of the expressions in Eq. ( j2.1| ) its full minimization with respect to p(x, z) for the 
boundary conditions described in Fig. [I] represents a big numerical challenge. Since we are 
mainly interested in the local interface position we refrain from seeking this full solution. 
Instead we restrict the space of possible density distributions to a subspace of piecewise 
constant densities (see Fig. |3|). Within this approximation at the position of the liquid- vapor 
interface z — l(x) the density varies steplike between the bulk values determined by Eq. fl2.8|) . 
Moreover, with p(x, z < d w ) = (see Fig. |3|) we take into account that there is an excluded 
volume near the substrate surface induced by the repulsive part of the substrate potential. 
Approximately one has d w = ~(o~f + o~ w ). Thus the analysis of the grand canonical functional 
amounts to inserting the steplike density distribution ( "sharp- kink approximation" ) 

p{x, z) = G{z - d w ){ Pl e(l(x) -z) + p g Q{z - l{x))} (2.10) 

into the functional ( |2.1| ), with l(\x\ — > oo) asymptotically approaching the function 

a(x) = l Q(-x) + (l + x tan 6)Q(x); (2.11) 

9 is the contact angle (see Fig. [TJ). 

With this ansatz the grand canonical free energy functional can be systematically de- 
composed into bulk, surface, and line contributions. The decomposition is carried out for 
a finite system, and in a second step the thermodynamic limit is performed. This lengthy 
calculation is carried out in detail in Ref. [|Tj]]. Here we quote only the main results. There 



are artificial surface and line contributions which stem from truncating the system at finite 
distances before considering the thermodynamic limit; we omit them here. One obtains the 
following expression for Q[p\: 

Q[p(x, z)\ = A«(0) Q b ( Pl ) + A® (6) n b ( Pg ) 

+Att s (l ,9)+L y n l [l(x)}, (2.12) 

with the volumes A w = \L x L y L(6) and A^ = \L x L y L(6), where 1(6) = L z - l = 
i-Lj; tan(0), and the surface area A = L x L y (see Fig. 0). The first two terms in Eq. ( |2.12| ) 
describe the bulk free energies of the liquid and vapor phases, with the bulk free energy 
density ^b(p-y) given by Eq. (|2.8| ). The surface contribution O s (Z , 0) consists of the following 
terms: 

Q S {1 , 6) = l Q b {pi) + o w i + \o ig + - 1 o ig + \uj(1 ) (2.13) 

2 ZCOS{o) Z 

where a w i and ai g are the wall-liquid and liquid-vapor surface tension, respectively. This 
expression corresponds to that for thin liquidlike wetting films adsorbed on homogeneous 
and planar substrates as obtained by the same approach ||. In this context it is shown that 
the equilibrium thickness Iq of the liquidlike film minimizes the effective interface potential 

( 7 7 \ 

= Ap | Pl J dzt(z) - J dzV(z) J = J2 J, 1 » ( 2 - 14 ) 
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a2 is known as the Hamaker constant. The interaction potential t(z) of a fluid particle with 
a half space occupied by fluid particles is given by 

oo 

t(z) = I dz' I d\ { w{Jr\ + z 12 ). (2.15) 

z 

The properties of the effective interface potential determine the character of the different 
wetting transitions. The contact angle 9 follows from Young's equation: 

a wg — o wl uj(Iq) , . 

cos 9 = — = 1 H (2.16) 

with the liquid-vapor surface tension given within sharp-kink approximation by 

oo 

a lg = -l(Ap) 2 [ dzt(z). (2.17) 



The line contribution can be split up into one term independent of l(x) and one func- 
tionally depending on l(x): 

n l [l(x)]=r a (l ,9) + r l [l(x)} (2.18) 

with 

r.M = ?(AP) 2 fl-AVA P. (2-19) 



4:71 \ ta.n9 / tan# 



2tt 



|[(/ - d w ) 2 + 07] (-^ + arctan(/ - d w )^j + l - (4 j + ^ j- ^ ~ + ' 

i>i Zl ' 



+ 1) 



and 

7?[/(a;)] =r u [Z(x)]+ri[Z(a;)] (2.20) 
where the first term is given by an integral over the effective interface potential: 

oo 

Tu[l{x)] = -Ap { pi J dx J dz (pit(z - d w ) - V(z)) 

-oo a (x) 



dx{u{l{x)) -u(a(x))} (2.21) 
with given by Eq. ( |2.14|) . The expression 
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tS(x)} = -(Ap) 



oo oo l(x) a(x') 

2 j dx J dx' J dz J dz' w(\x — x'\,\z — z'\) 

oo -oo a(x) -oo 



oo oo 



oo 

**J* J** l"*te-A,\*-A) 

oo -oo a (x) l(x') 



(2.22) 



describes the free energy contribution from the liquid- vapor interface and is a nonlocal func- 
tional of l(x). In Eq. (|2.22|) w is an integral over the attractive fluid-fluid interaction 



w( 



(x,z) = / dyw(\/x 2 + y 2 + z 2 ). 



^2.23) 



The expansion of the integrand in Eq. ( |2.22j ) into a Taylor series with respect to 



yields the gradient expansion of the functional expression, with the leading term 



^«*>i=<*/HHSW i+ (i)' 



(2.24) 



which is now a local functional of l(x). 

The line contribution (Eq. ( |2.18| )) to the free energy is minimized by using either 
Eq. (|2.22|) or Eq. (|2.24 ), respectively, yielding the equilibrium liquid- vapor interface l{x) 
and the corresponding line tension r as the minimum value fij [/(#)]. In the following the 
application of Eq. (|2.22|) will be called "nonlocal theory", whereas Eq. (|2.24|) is used for 
the so-called "local theory" , also known as "interface displacement model" . The expression 
in Eq. ( |2.24| ) measures the variation of surface free energy due to the deformation of the 
liquid-vapor interface. Equations (|2.18| )- (|2.21| ) and (|2.24j ) provide a prescription of how to 
introduce the parameters serving for a microscopic description of the underlying molecular 
interactions into the more simple interface displacement model which is originally motivated 
by phenomenological, macroscopic considerations. 



B. Intrinsic structure of the contact line 



First we analyze the shape of the interface within the local theory. Minimization of 
Eq. ( [2.18|) by using the functional in Eq. ( |2.24| ) leads to solving the corresponding Euler- 
Lagrange equation (ELE) following from Sfl[/Sl(x) = 0: 

= -Ap[ Pl t(l(x) - d w ) - V(l(x))] = ^'^ 2)3/2 = ^K(x). (2.25) 

K(x) is the curvature of the planar trajectory (x,l(x)), with its local radius of curvature 
- which is also one of the two principal radii of curvature of the surface (x,y, z = l{x)) - 
given by R = 1/K. It is related to the mean curvature H of the manifold (x,y,z = l(x)) 



dm (I) 
dl 
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by H(x) = ~K(x). Equation ( |2.25| ) relates the local curvature K(x) to the derivative 



of the effective potential (Eq. (|2.14f) ) exerted on the liquid- vapor interface. The equation 



is often referred to as "augmented Young-Laplace equation" [33]. This nonlinear ordinary 
differential equation can easily be solved using standard numerical tools, e.g., the algorithms 
from the Numerical Algorithm Group (NAG). We have used a Runge-Kutta algorithm with 
a prescribed starting point (x , Iq) and an initial derivative 1' . If 1' = 0, the solution is l{x) = 
lo = const, which describes the thin wetting film of constant thickness on the homogeneous 
and planar substrate. If, however, a very small initial derivative, e.g., 1' = 10~ 8 , is chosen, 
then the solution diverges for increasing values of the lateral coordinate x, asymptotically 
exhibiting the constant slope tan# for \x\ — > oo. The origin of the coordinate system is then 
shifted in lateral direction such that the position x = corresponds to the intersection of 
the asymptotes of l(x). 

Figure f| shows liquid- vapor interfaces around a three-phase contact line for a system un- 
dergoing a first-order wetting transition. Figure |(a) displays the temperature dependence 
of the profiles l(x), whereas Fig. |](b) presents the deviation of the profiles from the asymp- 
totes, i.e., 5l(x) = l(x) — a(x). Apart from a different choice of interaction parameters and 
thus a different wetting transition temperature, these results are in agreement with those 
obtained in Ref. [|17|. One of the main features in case of a first-order wetting transition is 
that l(x — > oo) approaches its asymptote from below (5l(x > 0) < 0), whereas l(x — > — oo) 
approaches its asymptote from above (8l(x < 0) > 0) (see also Ref. [p0|1 ). A more detailed 
discussion of the interface morphology is given in Ref. [[17]]. 

Within the nonlocal theory the ELE is the nonlocal integral equation 

oo l(x')-l(x) 

-Ap[p l t(l(x)-d w )-V(l(x))] = -(Ap) 2 J dx' J dz'w{\x-x'\,\z'\), (2.26) 

-oo 

i.e., the right-hand side of the local differential ELE ( |2.25| ) is replaced by an integral expres- 



sion. In Ref. [[17] the interface profiles within the nonlocal theory have been obtained from 
a numerical solution of the ELE (|2.26|) . The insertion of the discretization (xj, U = l(xj)) of 
the function l(x) on a lattice with N mesh points into Eq. ( p.26j ) yields a set of N coupled 
equations for the values U which can be solved using standard numerical tools. This in itself 
is already a demanding task compared to the relatively easy numerical solution of the local 
ELE ( p.25| ). But moreover the solution distinctly depends on the number of mesh points 
N, a fact that due to low computer power has not been revealed in Ref. |T7[. An enhance- 



ment of N such that the numerical resolution is satisfactory would require large computer 
memory and a very long time of computation. For this reason we numerically minimize the 
functional ilf 1 itself, i.e., Eq. ( j2.18| ) with the nonlocal expression (|2.22|) for T; [/(#)], 

instead of solving the ELE. To this end the interface is also discretized (xj, /j = l(xi)) with a 
linear interpolation between the mesh points. The z integration in Eq. ( [2.21 ) and the z and 
z' integrations in Eq. ( |2.22j ) can be carried out analytically. With the remaining integrations 
performed numerically, the discretization procedure yields a function Qi ({/«}) which is mini- 
mized with respect to the N variables U using a particulary suitable algorithm from the NAG 
routine library. The price of a significant increase of the numerical effort, which is required 
by this approach as opposed to the seemingly easier way of solving Eq. ( |2.26j ) directly, is 
justified by the fact that the solution does not depend sensitively on the choice of the mesh 
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size Sx = Xi + i—Xi, as it does for the ELE (|2.26|) . A repetition of the numerical minimization 
using a mesh with a larger number of points N yields an identical result, which corroborates 
the reliability of the minimization procedure. Moreover, our results are in accordance with 
those obtained for geometrically structured, wegde-shaped systems which also indicate that 
interface profiles obtained from the local and the nonlocal theory differ only slightly |[43| . 

The same minimization procedure can be applied to check the accuracy of the calculation 
of interface profiles within the local theory As expected, the minimization of the expression 
( [2.181 ) using the local functional Eq. ( |2.24| ) yields results identical to those obtained from 
the solution of the differential equation ( |2.25| ) whose numerical solutions do not depend on 
the number of mesh points. 

For the same system as discussed in Fig. |, Fig. || displays typical examples for the 
difference between the predictions obtained from the local and the nonlocal theory in terms 
of the deviation from the asymptotes 5l(x). Figure § shows the curvature K(x) of the 
respective interfaces. K(x) is slightly underestimated within the local theory. These figures 
clearly demonstrate that the obvious differences between the predictions of both theories 
are small and confined to those regions where K(x) is largest. Significant deviations of the 
nonlocal predictions from the local ones only occur at sections of the interface where the 
radius of curvature 1/K(x) is less than 10a f. This corrects the nonlocal results obtained in 
Ref. [17] from solving the ELE. Our refined numerical analysis shows that the differences 
are largest for high values of contact angle, i.e., for temperatures far below the wetting 
temperature. Upon approaching the wetting transition the contact angle and the curvature 
of the interface decrease, and the difference between the local and the nonlocal predictions 
becomes smaller. 

Figure ^ displays liquid-vapor interfaces calculated within the local theory for a system 
which exhibits a critical wetting transition at T* = 1.2. In the case of critical wetting l(x) 
approaches the asympotes from above both for x — > oo and for x — > — oo, i.e., 5l(x) is 
always positive. If the system undergoes a critical wetting transition, the contact angles are 
rather small, and the interface profile exhibits very small local curvatures (see Fig. |7[) and 
an extremely broad transition region (which is the interval where the deviation of l(x) from 
its asymptotes is larger than a certain prescribed value). 

In the case of a critical wetting transition it is presently impossible to perform a full 
numerical minimization of Q^ 10 ^: due to the huge required system sizes over which the 
numerical integrations in Eqs. ( [2.2 1| ) and (|2.22| ) have to be performed the available computer 
resources are by far exceeded. However, for a given equilibrium profile l(x) calculated within 
the local theory the relative difference (Clf 11 "^ —VLf oc ^)/Vtf oc \ with {l\ nloc ^ evaluated for l(x), 
is of the order of 10~ 3 . Therefore in the case of critical wetting we expect the difference 
between the interface profiles (as well as the line tensions) calculated within the local and 
the nonlocal theory to be negligibly small. 



C. Line tension 



The line tension r corresponding to the equilibrium liquid-vapor interface profile follows 
from Eqs. (glg)-(gg§: 



T «n)loc) = Q ((n)loc) ^ = Q ((n)/ c) ^ ^ 

{l(x)} 



9 



Within the local theory 

T (loc) is 

obtained by inserting the solution l(x) of the ELE into Q\ ; 
within the nonlocal theory r^ 100 ^ immediately follows from the minimization procedure itself. 

Figure |8| shows the results for r^ oc ) and r^ nloc ^ for the system undergoing a first-order 
wetting transition at T* = 1.102 (compare Figs. f|-|6|). Due to the large numerical effort only a 
few results have been obtained in the nonlocal theory. The data indicate that the difference 
T (nioc) _ T (ioc) - g rather sma ll and vanishes upon approaching the wetting transition. A 
thorough discussion of the different contributions to the line tension is presented in Ref. ||17||, 



including a comparison with the singular behavior predicted by Indekeu |[19| , |20[| . The results 
obtained in Ref. flD]] for the local theory are still valid; therefore we can refrain from a further 
discussion here. We only show in Fig. || the temperature behavior of the two contributions 
and Ti (Eqs. ( |2.20[ )-( p.24 )) to the line tension r, as derived within local and nonlocal 
theory. The difference jf 110 ^ — rf oc ^ is about twice the difference t^ 1,1 "^ — t^ oc \ and both 
vanish for T — > T w . 



III. ADSORPTION ON A CHEMICAL STEP 
A. Model and density functional theory 

The study of the system shown in Fig. |] is a first step towards the understanding of 



wetting on chemically heterogeneous surfaces P7| . In this case the translational symmetry 



in lateral directions is broken by the variation of the substrate potential in x direction 
(Fig. 0). We describe the fluid in the half space {r G R 3 \z > 0} as in Sec. (Eq. Q). The 
only difference is that the substrate potential V(z) in Eq. ( [2.2| ) is replaced by V(x, z). As 
boundary condition one has p(x, z — > oo) = p g for all x, including \x\ — > oo. To be specific 
we assume that the substrate potential is the pairwise superposition of Lennard- Jones pair 
potentials between the fluid and substrate particles, 



Ur) = 4e : 



<t±\ 12 _ (^±y 



J \ r J 



(3.1) 



where the + and — signs stand for interaction of a fluid particle with the different species 
occupying the quarter spaces w + = {r £ IR 3 |x > 0, z < 0} and if_ = {r G R 3 |x < 0, z < 0}, 
respectively. This superposition yields [^7[] 

u f + u 3 1 u t ~ u z ( 1 ( r \ 3 3 1 



Vatt(x,z) = -^-^^ + ^—^ -- - + 



2 z 3 2 \x 3 \xzJ Ixzr 

U^ z + V% z 1 U^ z - Ul z fx 1 X 



+ t:- 



2 z 4 2 \z A r 2z 2 r 3 

+^^^fA-f4 + ^V^- m ^ rn + n>5) (3.2) 
2 V x 4 V x 4 r 2 x^r 6 I I 



for the attractive part of the substrate potential, where r = r(x, z) = \/ x 2 + z 2 . The 
coefficients uf, uf z , and functions of the interaction potential parameters e± and 

a±; uf and uf z are the coefficients of the corresponding homogeneous, flat, semi- infinite 
substrate w± occupied by the species "±" with the substrate potential of a homogeneous 
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substrate defined as in Eq. (|2.2|). For the repulsive contribution to the substrate potential 
we use the simple ansatz 



v rep (x, z) = e(-x)=l + e(x)4, (3.3) 

Z J Z J 



because the repulsive interaction of each half of the substrate decays rapidly. 

This system has been studied in detail in Ref. |27|. As stated in the Introduction we 
revisit this problem in order to study the interface morphology with the refined numerical 
analysis described in Sec. |J. Moreover, in Ref. [[27]] only the cases of critical and complete 
wetting have been covered. Here we also study a first-order wetting transition for which one 
can expect the largest differences between the local and the nonlocal approach (see Sec. 0). 

The behavior of the substrate potential V(x,z) for large \x\ and z fixed, i.e., far from 
the heterogeneity, 



V (x - ±00, z) = -4 - 4 + 4 ± + °(W" 4 )' ( 3 - 4 ) 

z 6 z 4 z y \xr 



causes the liquid- vapor interface to asymptotically approach the constant values l(x — > 
±00) = l± valid for the corresponding homogeneous substrate. Therefore it is appropriate 



to make the following sharp-kink ansatz for the particle density distribution (see Eq. (|2.10|) ): 



p(x, z) = [Q(-x)Q(z - d-) + Q(x)Q(z - d+)} 

x [ Pl Q(l(x) -z) + Pg Q(z - l(x))\. (3.5) 

For this density distribution the grand canonical free energy functional decomposes into 
bulk, surface and line contributions such that the surface contribution describes the wetting 
of each single homogeneous, flat, semi-infinite substrate: 

n([p(x,z)};T,p; [w], [V]) = A{l b (p g , T, p) + A£l s (l±; T, p; [w], [V]) 

+L y n l ([l(x));T,p:-[w},[V}), (3.6) 



where fib is given by Eq. ( |2.8| ), A = L x L y L z is the volume filled with fluid, and A = L x L y is 
the surface area of the substrate surface (see Fig. ||). Again artificial terms due to truncations 
are omitted. 

to,(l±) = -(/+ + l-)AQ b + dig + -(a wl (l+) + a w i(l-)) + -(w(Z+) + u(Q) (3.7) 

is the arithmetic mean of the surface free energy densities corresponding to wetting of the 
homogeneous substrates w + and W- by liquidlike layers of thickness l + and /_, respectively. 
AQb = ApA/i + 0((Ap) 2 ) with Ap = pi — p g measures the undersaturation if the liquid 
phase is not yet thermodynamically stable, i.e., if Ap = po — p > 0. The line contribution 
fi;[/(x)] reads 

Qi[l(x)] = T wl + r lg + u[l(x)] (3.8) 

with 
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Zoo (a:) 



oo oo loo{x) — di, CK 

*[<(*)] = A* / «fc(I(,) - Mx» + APP, / «fc / *«« - Ap/ «fe / 

-°° -°° l{x)-d+ -oo 



i(a0 



/ 



c/.r 



dzt(x, z) 



-oo 



I (a:) — 

oo oo oo Z(jc)— Z(ai') 



oo Zoo(^)— d« 

da; 

x l(x)-du 



dzt(x, z] 



I 



— oo — oo 

where Zoo(^) = Q(—x)l- + 0(x)/ + and 



oo oo 

' / Jj,-7.f„l J\ 



t(x,z) = / dx / dz w(x', z). 



(3.9) 



The rather lengthy expressions T wi and T; 5 do not depend on the function \{x) and are 
given in Ref. P7| ]. w[i(a;)] is a nonlocal functional which has to be minimized to yield the 
equilibrium profile l(x) and the line tension r = f^[Z(a;)]. 

The local approximation of the quadruple integral in u)[/(a;)] - which describes the free 
energy contribution due to deformation of the liquid-vapor interface - is given by the local 
functional 



oo oo oo 



l(x)-l(x') 



-far J to J "J* J <**Q*-*\,\'-A 



-oo — oo 
oo 



°ic \ dx< \ 1 



dl 

dx 



(3.10) 



B. Numerical results for the interface morphology 

Within the local theory the interface morphology follows from the ELE 

Atl b - Ap[ Pl t(l(x) - d w ) - V(x,l(x))} = f^M^f, = *i 9 K(x). (3.11) 

Apart from the finite undersaturation AQ^ > and the different substrate potential the 
structure of this ELE is the same as in Eq. ( |2.25| ). For reasons of simplicity we assume that 
d w = d+ = d~. Guided by the experience from Sec. |I| the equilibrium profile within the 
nonlocal theory is obtained by minimizing Qi (Eq. ( |3.8| )) instead of solving the corresponding 



ELE as it has been done in Ref. [p^] . Here we apply the same numerical procedures as in 
Sec. 0. 
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Figure |]Xj displays interface profiles on a complete wetting path, i.e., along an isotherm 
(T = const, A/x — > 0); the system exhibits critical wetting transitions at = ksT w /tf = 1.2 
on the substrate w_ and at T* = 1.0 on the substrate w + , respectively. For the choice 
T* = 1.1 at coexistence the substrate w + is completely wet whereas u>_ is only partially 
wet. We find that the difference between the local and the nonlocal results for l(x) as well as 
for the line tension r is smaller than the numerical resolution and therefore not significant. 
The minimization procedure requires a large numerical effort and is only applied for values 
of \i not too close to coexistence /i = /x such that the effective size of the system, which has 
to be treated numerically, is not too large. 

For larger effective system sizes, i.e., close to coexistence A/i = 0, we have to resort 
to an alternative procedure for evaluating the equilibrium interface profile l(x) within the 
nonlocal theory. The corresponding ELE 

oo l(x')-l{x) 

AVt b - Ap[p l t(l{x)-d w )-V{x,l{x))} = -{Apf J dx J dz'w(\x-x'\,\z'\) (3.12) 

-oo 

(compare Eq. ( |2.26| )) is solved numerically for a certain choice of thermodynamical pa- 
rameters (T,/i). To this end the lateral extension of the system is truncated and l(x) is 
discretized on a lattice with iV lattice points and a mesh size Sx, yielding a set of N coupled 
equations for the values U = l(xi). These equations are solved using a suitable numerical 
algorithm. The solutions depend on the mesh size 5x in such a way that solutions l(x; 5x) 
for different 5x differ only with respect to a simple rescaling of the x axis by a constant 
numerical factor. We note that this particular scaling behavior of the numerical solutions of 
the ELE only occurs in the present case of a substrate with a chemical step but not for the 
intrinsic structure of the contact lines investigated in Sec. [TJ]. Thus in the latter case this 
alternative approach is not applicable. The procedure of finding the correct solution is as 
follows. We calculate solutions of the ELE on lattices with different Sx. A suitable length 
scale is defined which measures the typical width A on which the profile varies between its 
asymptotic values l + and /_, e.g., the distance between the points where l(x) deviates by 
10% from the asymptotes. As a function of 5x the width A asymptotically approaches a 
linear dependence A(5x — > 0) = A 5x + B which can be extrapolated to 5x = 0. This yields 
a factor 70 = A(5x = 0)/A(&r ) by which a certain solution l(x; Sxq) of the ELE has to be 
scaled in x direction in order to find the limiting profile corresponding to an infinitely fine 
lattice (see Fig. [TI]). In order to test the reliability of this approach we have carried out the 
following additional cross-check: a solution l(x; 5xq) of the ELE obtained by using a mesh 
size Sxq is scaled in x direction by different factors, i.e., x 1— > 71; the correct scaling factor 
70 is the one for which the line contribution to the free energy f2/[/(7x; Sxq)] is minimized. 
Both procedures are based on the assumption that the true solution differs from each of the 
numerically obtained solutions l(x; 8x) only by a simple rescaling of the x axis. (For those 
systems, for which the full minimization and the solution of the ELE can be carried out 
both, this assumption has been verified.) Within the numerical accuracy the two methods 
give identical results for the scaling factor. Moreover, as expected the solution obtained by 
rescaling of the x axis with this optimal factor is indistinguishable from the one predicted 
by the local theory. 

Figure shows interface profiles for a thermodynamic path along coexistence (T — > 
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T w ,A/i = 0). The system exhibits critical wetting at the same transition temperature 
T* = 1.0 for both substrates w + and the parameters of the substrate potential are chosen 
such that the asymptotic equilibrium film thicknesses l + and /_ are different although the 
wetting transition temperatures are equal. Due to the broad crossover region these profiles 
have been obtained by using the scaling procedure as described for finding the minimum 
of Qi. Also in this case the local and the nonlocal results are indistinguishable. Not only 
the profiles but also the line tensions obtained within the local and the nonlocal theory, 
respectively, are very close. The relative difference of the line tensions (r^ n ' oc ^ — r^ oc ')/T^ oc ^ 
is of the order of 10~ 4 , both for complete as well as for critical wetting. 

In Fig. ^ we present results for interface profiles for a system with very thin liquidlike 
wetting films. Both corresponding homogeneous substrates undergo first-order wetting tran- 
sitions: w + at the temperature T* rs 1.102 and W- at m 1.314; both wetting temperatures 
are above the temperatures considered here. In the previous examples the differences be- 
tween the results obtained within the local and nonlocal theory are small. But here they are 
detectable since the interface curvatures are larger than in a system with thick wetting films 
as they occur for critical and complete wetting. The relative difference (r^ n ' oc ^ — r^ loc ') / t^ 1oc ' 
is of the order of 10 -2 . 



IV. SUMMARY 

We have analyzed the morphology and the associated line tensions of liquidlike wetting 
films at a three-phase contact line (Fig. [I]) and on a planar substrate across a chemical 
step (Fig. |2|). By using refined numerical techniques we have compared quantitatively the 
predictions obtained within a local displacement model for the interface profile and within a 
nonlocal density functional theory. Based on general arguments p3|-[26[] the latter approach 



is the more accurate one. We have obtained the following main results: 

1. Within the present mean-field theories the equilibrium interfacial profiles are deter- 
mined as the minimum of the line contribution to the functional of the grand canonical 
free energy. These equilibrium profiles, within the sharp-kink approximation (Fig. [|), 
can be obtained either by an explicit minimization procedure or by solving the cor- 
responding Euler-Lagrange equation for the minimum. Within the local theory both 
approaches are equally successful and robust. However, within the nonlocal theory in 
practical terms only the minimization procedure yields reliable access to the true mini- 
mum whereas solving the necessarily discretized version of the nonlocal Euler-Lagrange 
equation leaves one with a very slow convergence with respect to the mesh size of the 
discretization. This mesh size dependence of the solutions within the nonlocal theory 
was not revealed in Refs. |H] and |27j and led to a significant overestimation of the 
quantitative discrepancy between the predictions of the local and the nonlocal theory. 

2. Concerning the intrinsic structure of the three-phase contact line the local and the 
nonlocal theory yield only slightly different interface morphologies (Figs. |] and |^) and 
line tensions (Figs. [8] and |^) in the case of first-order wetting. In the case of critical 
wetting the two theories yield indistinguishable results (Fig. [7]). For first-order wetting 
the deviations between the results obtained from the local and the nonlocal theories are 
confined to regions within which the radius of curvature of the interface profile is less 
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than about 10 atomic diameters 07 of the fluid particles (Fig. ^|). Upon approaching a 
first-order wetting transition the line tension r of the three-phase contact line diverges 
logarithmically (Figs. |8] and |9|). r changes sign as a function of temperature and 
typically is of the order of e//er/ where ej is the depth of the interaction potential 
between the fluid particles. 

Concerning the morphology of a wetting film on a planar substrate with a chemical 
step the results obtained from the local and the nonlocal theory for complete and 



critical wetting are indistinguishable (Figs. |T0| and |12|) . Small differences appear only 
in the case of first-order wetting (Fig. 



Thus we conclude that the local interface displacement model yields quantitatively reliable 
results for the interface morphology and the line tension of laterally inhomogeneous wetting 
films. The more accurate nonlocal description is necessary only for cases in which the 



interface profiles exhibit large curvatures. Thus all results and discussions in Refs. |L7 
and |27[ based on the local theory turn out to provide actually a quantitatively reliable 
description of the systems under consideration. 
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FIGURES 




FIG. 1. Shape l{x) of the liquid-vapor interface near the three-phase contact line on a homo- 
geneous, planar substrate whose surface is located at z = 0. is the thickness of the microscopic 
liquid wetting film, d w denotes an excluded volume due to the repulsive part of the substrate po- 
tential. a(x) is the asymptote l(x) in the limit |x| — > oo. The angle 9 between a(x > 0) and the 
substrate surface is the contact angle given by Young's law cos# = (a wg — a w i)/ai g . The inter- 
section between the asymptotes a(x < 0) and a(x > 0) defines the position x = 0. In order to 
facilitate the proper thermodynamic limit the system is truncated at x = ±L x /2 and z = L z . The 
configuration is taken to be translationally invariant in the y direction. 
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FIG. 2. Morphology l{x) of a wetting film which covers a planar substrate with surface z = 0. 
The substrate consists of two halves meeting at x = 0. The material filling the left (right) half 
favors a thin (thick) wetting film, = l{x — > =Foo) are the equilibrium film thicknesses of the 
corresponding homogeneous substrates, which are characterized also by different excluded volumes 
dj. The system is translationally invariant in the y direction. It is truncated at z = L z and 
x = ^L x /2 in order to facilitate the proper thermodynamic limit. 
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z=d 



w 



FIG. 3. Sharp-kink approximation for the particle density distribution. At the position 
z = l{x) of the liquid- vapor interface the particle density distribution varies steplike between the 
constant liquid (pi) and the constant gas density (p g ). At z = L z the density is truncated in 
order to facilitate the thermodynamic limit. Moreover, the density vanishes for z < d w due to the 
repulsion between the fluid and the substrate particles. 
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FIG. 4. Temperature dependence of the shape of the liquid-vapor interface calcu- 
lated within the local theory for a system which exhibits a first-order wetting transition at 
T* = kBT w /ef ps 1.102. The substrate potential parameters are d w = 1.05cj/, 113 = 3.710e/<7j, 
U4 = 5.566e/Oj, and ug = 0.876e/Oj. In (a) the interface profiles l(x) (full lines) and their asymp- 
totes a(x) (dashed lines) are shown, whereas (b) displays the deviation 51 (x) = l{x) — a(x) of l(x) 
from its asymptotes. The corresponding contact angles are 6 ps 44.5°, ps 36.9°, and 6 ps 26.6° for 
T* = 0.8, 0.9, and 1.0, respectively. l(x) approaches its asymptotes from above for x — ► — 00 and 
from below for x — > 00. l(x) closely follows the asymptotes for large values of the contact angle, 
i.e., for low temperatures. 



21 




x/a f 

FIG. 5. Comparison between liquid-vapor interface profiles calculated within the local (dashed 
lines) and the nonlocal (full lines) theory for the same system as discussed in Fig. f|, in terms of 
the difference 5l(x) = l(x) — a(x) for T* = 0.8 (a) and T* = 1.0 (b). The difference between the 
local and the nonlocal results decreases upon approaching the wetting transition. The local theory 
slightly underestimates the local curvature K(x). 
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FIG. 6. Comparison between the curvatures K{x) = I "(x) /(1 + (I '(x)) 2 ) 3 / 2 of the liquid- vapor 
interface profiles as calculated within the local (dashed lines) and the nonlocal (full lines) theory for 
the same system as in Fig. ||. The parameters chosen for (a) and (b) are the same as for Figs. ||(a) 
and ^(b), respectively. Also for the curvature the difference between the local and the nonlocal 
results decreases upon approaching the wetting transition. The local theory slightly underestimates 
the curvature. 
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FIG. 7. Liquid-vapor interfaces calculated within the local theory for a system which exhibits 
a critical wetting transition at T£ = 1.2. Here the substrate potential parameters are d w = cry, 
us = 1.683eyCTp = 10.098e/Oj, and ug = 0.224e/Oj. In (a) the interface profiles l(x) (full 
lines) and their asymptotes a{x) (dashed lines) are shown, whereas (b) displays the deviation 
5l(x) = l(x) — a{x) of l(x) from its asymptotes. The corresponding contact angles are 9 ~ 0.31°, 
9 rj 0.17°, and 9 m 0.06° for T* = 1.14, 1.16, and 1.18, respectively. l(x) approaches its asymptotes 
from above both for x — > —oo and x — > oo. In (c) the curvatures K(x) are plotted. As compared 
to the system exhibiting a first-order wetting transition (see Figs. [| - ||) 5l(x) and the radii of 
curvature 1/K(x) are extremely large . Due to numerical difficulties no data for the nonlocal 
theory are available for this system. However, since the interface curvature is several orders of 
magnitude smaller than in the system which displays a first-order wetting transition - with already 
small deviations of the local from the nonlocal theory - the difference between local and nonlocal 
results for the present system is expected to be indistinguishable on this scale. 
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FIG. 8. Line tensions as predicted by the local and the nonlocal theory for the same system 
as in Fig. ||. The transition temperature for first-order wetting is ~ 1.102. The difference 
T (nloc) _ T (loc) vanishgg U p Q n approaching T w . The singular behavior of the line tension can be 
derived analytically within the local theory; upon approaching the wetting transition r diverges as 
t(T — ► T w ) ~ ln(0(T)), as indicated in the inset. (The dotted line is the estimated extrapolation 
according to the logarithmic divergence.) We note that as function of temperature r changes sign. 
Apart from the close vicinity of T w and its zero the absolute value of r is of the order of e//c/- 
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FIG. 9. Temperature dependence of the contributions r w and t% (Eqs. ( |2.20| )-( |2T2~4| )) to the 
line tensions r( n ' oc ) and r(' oc ) for the same system as in Figs. || and ||. For T — > T w ~ 1.102ej/fes 



both t. 



(nloc) 



-T { j oc) andr (nioc) 



_(Zoc) 



r ' vanish. As already discussed in Ref. [17], the leading singular 



behaviors t u ~ —1/6 and Tj ~ 1/0 for 8 —* 0, i.e. 
leaving the residual singularity r ~ — ln(0). 



T —>■ T w compensate each other in r = r w + Tj, 
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FIG. 10. Shape of the liquid-vapor interface on a planar substrate exhibiting a chemical step 
(see Fig. ||) within nonlocal and local theory along an isotherm (T* = 1.1, A/u* = Af^/ej — ► 0), 
i.e., for complete wetting. The parameters for the substrate potential are chosen such that both 
substrates individually exhibit critical wetting transitions, the substrate W- (x < 0) at = 1.2 



and the substrate w + (x > 0) at T£ = 1.0: 



d+ 



I In 



2.079e/cr^, 



1.2 



and Ug = 0.277e/O"j, whereas the parameters d w and u^- are chosen as in Fig. |7] with ti 4 ^ it 



12A75e f aj, 

A,z- 

Therefore for this temperature at coexistence A/i = the substrate w+ is wet whereas the substrate 
W- remains only partially wet so that l(x — > oo, A/i = 0) = oo and l(x — * — oo, Afj, = 0) < oo. The 
nonlocal result for A[i = 1.5 • 10 -6 has been obtained by solving the discretized version of the ELE 
on lattices of different mesh size 5x and the application of a pointwise extrapolation scheme 5x^0 
as outlined in the main text. The other interfacial shapes calculated within the nonlocal theory 
result from a free numerical minimization of the line contribution to the free energy 
The differences between the local and the nonlocal results are not visible on this scale. Therefore 
only one type of line is shown for each interfacial profile. 
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FIG. 11. Dependence of the numerical solutions of the nonlocal ELE on the mesh size 5x 
for the same set of parameters as in Fig. 10 and for Afi* = 10~ 3 . (a) displays profiles obtained 
for different values of 5x. As 5x is decreased, the distance A between the two points where l(x) 
deviates from the asymptotes /+ and /_ by 10% becomes smaller, as shown in (b). This dependence 
can be extrapolated to the limit 5x = 0, as indicated by the straight line in (b). Using the value of 
A from this extrapolation, the correct result in the limit of an infinitely fine lattice can be obtained 
along the lines discussed in the main text. This limiting profile is shown in (a) {5x = 0); it turns 
out to be indistinguishable from both the local result and the result obtained by the full numerical 



minimization of Q, 



(nloc) 
I 



which can also be carried out for the present example. 
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FIG. 12. Liquid-vapor interfaces on a planar substrate with a chemical step (see Fig. |2|) as 
obtained within the nonlocal and the local theory. The parameters for the substrate potential 
are chosen such that both substrates exhibit critical wetting at the same transition temperature 
T* = 1.0. The parameters d+ and u+ are chosen as in Fig. 10, whereas d~ = aj, u% = u^, 
Ug = Ug, and u\ z = u\ x = 15.594e/ <r*. The profiles are calculated for different temperatures T* 
on a thermodynamic path along coexistence A[i = 0. The nonlocal results for both temperatures 
have been obtained by applying the extrapolation scheme explained in the main text. On the 
present scale the local and nonlocal results are indistinguishable and therefore only one type of 
lines is shown. 
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FIG. 13. Liquid-vapor interfaces on a substrate with a chemical step at x = (see Fig. as 
obtained within the nonlocal and the local theory. The parameters for the substrate potential are 
chosen such that both substrates individually exhibit first-order wetting transitions, the substrate 
W- at 1.314 (with d w = aj, = 2.513e/<7j, u^ z = u^ x = 3.770e/<Tp and Ug = 0.335e/crp 

and the substrate w + at « 1.102 (with the parameters d+ and u- as for Fig. ||] and u\ x = u\ z ). 
The system is at coexistence A/i = 0. Within the nonlocal and the local theory both profiles are 
obtained by a full numerical minimization of the line contribution to the free energy 
differences between the local and the nonlocal results are very small but detectable. The relative 
difference in line tensions (r^ ^ — t^ oc ))/t^ oc ) is of the order of 1CT 2 . 
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